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ABSTRACT 


Presented  in  this  report  are  the  results  of  an  investigation  to 
determine  analytically  two  thermodynamic  transport  properties  -- 
difference  in  specific  heats  and  the  change  in  entropy  --  as  described 
by  five  independent  equations  of  state  for  gases. 

Change  in  entropy  and  specific  heats  of  gases  are  used  in  thermo¬ 
dynamics  work  to  evaluate  properties  of  gases  undergoing  change  of 
state  such  as  flow  properties  behind  shocks.  Change  of  entropy  is  a 
logical  starting  point  for  this  work  and  where  applicable  properties 
are  known  at  high  temperatures  these  equations  may  be  used. 

Steam  was  used  for  illustrative  purposes  because  it  is  a  well 
defined,  thoroughly  explored  gas. 

The  development  of  the  equations  precedes  each  equation  in  a 
logical,  straightforward  manner  and  ends  with  the  desired  properties 
in  terms  of  known  thermodynamic  coordinates. 

Comparison  of  the  calculated  change  in  entropy  is  presented  in 
bar-graph  form  in  Figures  1  and  2  for  these  equations  of  state  for  a 
real  gas,  superheated  steam.  A  comparison  of  the  difference  in  specific 
heats  for  these  equations  for  the  same  gas  is  given  in  Figure  3. 

In  both  the  entropy  change  and  the  difference  in  specific  heat, 
good  agreement  is  found  among  the  five  equations  of  state. 
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I.  INTRODUCTION 

The  purpose  of  this  report  is  to  determine  the  dif¬ 
ference  in  specific  heats  at  constant  pressure  and  constant 
volume,  the  ratio  of  these  two,  and  the  change  in  entropy 
for  five  equations  of  state.  These  are  the  perfect  gas, 
Callendar,  Dieterici,  Berthelot,  and  Van  der  Waals  equations. 
After  these  relations  have  been  established,  some  examples 
of  how  they  might  be  used  will  be  shown  for  superheated 
steam. 

II.  DISCUSSION 

It  is  desirable  to  obtain  the  difference  in  specific 
heat  at  constant  pressure  and  the  specific  heat  at  constant 
volume  for  a  substance  in  terms  of  pressure,  volume,  and 
temperature. 

Heat  capacity  may  be  defined  as  the  ratio  of  heat 
.entering  the  system  and  the  change  in  temperature  of  the 
system  or  substance.  The  instantaneous  heat  capacity,  or 
specific  heat,  is  defined  as  this  ratio  as  both  the  heat 
added  and  temperature  change  become  smaller  and  smaller. 
Written  in  mathematical  terms,  _  ^  .  This  has  a  unique 

dT 

value  when  the  pressure  or  volume  of  a  substance  is  kept 
constant;  thus,  /dQ  \  and  (^\  *  This  requires, 

"  ^dTpj  ■  ^dTvj 

then,  that  specific  heat  be  a  function  of  at  least  two 
variables  in  thermodynamics. 

For  two  reversible  paths  joining  two  equilibrium 
states  of  a  system  from  Clausius’  theorem  we  may  write: 


Rj  Rj 


where  Rj^  and  R2  refer  to  the  path  taken  from  and  to  points 
i  and  f  in  any  thermodynamic  system. 
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which  means  that 


is  independent  of  the  path  taken 


and  allows  a  function  to  exist  in  the  coordinates  of  the 
system  whose  final  value  minus  the  value  at  the  initial 
noin-l-R  pniisl  s  the  intearel  .  Ue  mav  nail  it  the 


to  note  that  although  dQ  is  not  an  exact  differential,  ds 
is  exact.  The  reciprocal  of  the  Kelvin  temperature  is  the 
integrating  factor.  Dividing  both  sides  by  dT  and  multi¬ 
plying  by  T  gives  .  If  we  restrict  this  expression 

Hr  =  c!T" 

to  constant  volume  and  constant  pressure,  respectively,  we 
get; 


To  show  entropy  as  a  function  of  two  thermodynamic 
variables,  we  will  use  the  following  two  theorems  in  mathe¬ 
matics. 


Theorem  1.  If  a  relation  exists  among  x,  y,  and  z ,  we 
may  imagine  z  expressed  as  a  function  of  x  and  y;  whence, 

“  (ll)y 
If  we  let; 

^  ’  Ox)y  ■  ("^Ix  *  t^en,  dz  =  11  dx  +  IJ  dy, 

where  z,  M,  and  II  are  all  functions  of  x  and  y.  Differ¬ 
entiating  H  partially,  with  respect  to  y,  and  N  with  respect 
to  X,  we  get; 

/ aM \  _  a^z 

v'Sy/  X  ~  a  X  ay 

/aN^  .  3^z 
\7x/  y  *  ay  ax  * 


Since  the  two  second  derivatives 
it  follows  that: 


on 


(s) .  ■  (fi) 


the  right  are  equal. 
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This  is  known  as  the  condition  for  an  exact  differential. 


Theorem  2 .  If  quantity  f  is  a  function  of  x,  y,  and  z» 

and  a  relation  exists  among  x,  y,  and  z;  then  f  may  be 

regarded  as  a  function  of  any  two  of  x,  y,  and  z.  Similar¬ 
ly,  any  one  of  x,  y,  and  z  may  be  considered  to  be  a 

function  of  f  and  one  other  of  x,  y,  and  z.  Thus,  regard¬ 

ing  X  to  be  a  function  of  f  and  y, 

(t)y  (!7)f 

Considering  y  to  be  a  function  of  f  and  z, 

=  (lf)z  *  (lf)f 

Substituting  this  expression  for  dy  in  the  preceding 
equation,  we  get: 

"""•  ‘'’‘=  (lr)z  (#)f 

Equating  the  dz  terns  of  the  last  two  equations,  v.’e  get: 

It  may  be  shown  that  / dQ\  “  Cv  and  / dQ\  =  Cp. 

\dryv  \dT/P 

To  arrive  at  the  difference  in  specific  heats  in  the 
pressure,  volume,  and  temperature  coordinates  V7e  may  use 
entropy  as  a  function  of  T  and  V,  and  T  anti  P,  respective¬ 
ly. 

ds  =  dT  ♦  'lultiplying  through  by 

yTT^v  \7vyT 

T,  Tds  =  T  ^1^^  dT  ♦  dv  s  Cv  dT  t  T 
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From  the  equation  of  the  Helmholtz  function  A  =  U  -  Ts, 
we  get  the  differential  resulting  in  an  infitesimal  change. 
dA  =  dU  -  Tds  -  sdT  =  -pdv  -  sdT,  this  is  an  exact  differ¬ 
ential  since  A  is  an  actual  function.  From  the  above  mathe¬ 
matical  theorems  we  see  that  -/ap\  (ds\  Now: 

\^3Tjv  ■  y5vjT 


Tds  =  Cv  dT  +  T 


\3T/^ 


dv 


T. 


Imagining  entropy  to  be  a  function  of  T  and  p  v/e 

=  (w)p  *  (I|)t 


obtain:  _  /as\  .  /3s\  Multiplying  through  by 


Tds  =  T 


/  3s\  dT  +  T  (l£\  dp  =  Cp  dT  +  T 
\3T/p  VPr  V^)' 


dp. 


From  the  change  in  Gibbs  function,  G  =  H  -  Ts ,  dG  =  dH  - 
Tds  -  sdT  dG  =  vdp  -  sdT,  an  exact  differential.  There¬ 
fore,  /  \ 

(^)p  =  “(^)t  =  Cp  dT  4  T  (  -|v)  dp 


3T  p 

(  3V) 


Tds  =  Cp  dT  -  T  p  dp. 

Equating  these  two  equations:  Tds  =  Cp  dT  -  T  dp 

3T  p 

Tds  =  Cv  dT  +  T  dv 

Cp  dT  -  T  dp  2  Cv  dT  +  T  (3_p)  dv  i  solving  for 

3T  p  'aT  V 


dT,  CpdT  -  CvdT  =  T  ^|y^p  dp  +  T 


dv 


dT  = 


T  ( 3v) 


I  (^)  Imagining  T  to 


C^v  sT  p  ^  Cp-Cv  aT 
be  a  function  of  V  and  p,  dT  =  ^"^^p 


so 


CjT)  ^  (^)  -  T  (3V)  dp  +  T  (3p 

av  p  3p  V  Cp-Cv  IT  p  Cp-Cv  IT 


^  dv, 

V 
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Therefore,  the  terms  corresponding  to  dv  and  dp  on  each 
side  of  the  equation  are  equal. 


(3T)  _  T  (32)  (3T)  _  T  (_^) 

Tv  p  ■  Cp-Cv  Tt  V  3P  V  ■  Cp-Cv  3T  p 

Cp-Cv  =  T  ^12.5  and  Cp-Cv  =  T 

3T  V  3T  p  3T  p  IT  V 

which^  of  course,  is  the  same  from  each  equation. 

To  show  clearly  that  Cp-Cv  is  always  positive  for  all 
known  substances  let  us  get  the  change  in  pressure  with 
respect  to  volume  instead  of  (^p)  in  this  equation. 

Tt  V 


To  do  this  let  us  use  the  two  mathematical  theorems 
discussed  earlier.  Suppose  a  relationship  exists  between 
pressure,  P,  volume,  V,  and  temperature,  T,  so  that 
f(P,  V,  T)  =  0.  Let  us  imagine  P  as  a  function  of  v  and  T. 


dp  = 

(32) 

Tv  T 

dv  + 

(^) 
3T  V 

dT  and  v  as  a  -function  of  p  and  T. 

dv  = 

(3v) 

dv  + 

(32) 

dT  and  v  as  a  function  of  p  and  T. 

3P  T 

3T  V 

dv  = 

( av) 
ap  T 

dp  + 

(  3v) 
3T  p 

dT 

dp  s 

(32) 

Tv  T 

» 

(  3V )  (jp  • 

3p  T 

*■  (av)  dT 

TT  p 

• 

+  (^)  dT 

TT  V 

^41.5  dp  +  dT  +  dT 

3v  T  Tp  T  ^  3v  T  3T  p  3T  v 


dp  = 


Of* 


(i£) 


(av) 


_  +  (^) 

3v  T  TT  p  TT 


dT. 


Now,  of  the  three  coordinates  only  two  arc  independent. 
Let  us  choose  p  and  T  as  the  independent  coordinates.  The 
above  equation  must  be  true  for  all  values  of  dp  and  dT. 
First,  let  us  make  dT  =  0,  and  dp  not  equal  0. 


Then,  dp 


(3p)  (av) 
Iv  T  ap  T’ 


dp  ♦  0 
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If  dp  =  0  and  dT  not  equal  0  then, 


0 


( 9p)  ( 3v)  0  + 

av  T  3p  T 


(3v)  +  ( 3p) 

3V  T  IT  p  bT  V 


dT, 


so  then,  the  coefficient  of  dT  must  equal  0. 


( 3p)  (3v)  ( 3p) 

3v  T  TT  p  '3T  V 


0 


(3p)  (3v)  _  -(^) 

3v  T  TT  p  “  sT  V 


the 


V  “  Tv 
specific 


T  p  *  substituting  this  relation 

heat  difference  equation  we  get: 


into 


Cp-Cv  =  T  =  T 

Tr  V  Tr  p 


Cp-Cv  =  -T 

3T  p  3vT 


(^)  (3v)  (3V) 

“Sv  T  TT  p  P 


First  of  all, 


it  may 


be  remembered  that 


( 3p) 

Tv  T 


is 


always 


negative  for  any  known  substance  so  the  two  negative  signs 
when  multiplied  will  give  a  positive,  T  will  always  be 

positive  and  is  squared,  and,  whether  positive  or 

3T  p 

negative,  the  square  will  be  positive.  Therefore,  Cp-Cv 
will  always  be  positive,  or  Cv  is  never  greater  than  Cp, 


It  may  be  noted  that  as  T  approaches  absolute  zero  the 

( 3v) 

heat  capacities  are  equal,  Cp  may  equal  Cv  when  TT  p  =  ° 

as  in  the  case  of  water  at  maximum  density.  Next,  we  may 
determine  the  ratio  of  specific  heat  as  a  function  of 
P,  T,  and  V, 

From  the  two  Tds  equations  derived  previously. 


Tds  s  CpdT 
at  constant 


.  T  dp  ;  Tds  =  CvdT  ♦  T  dv 

P  3T  V 

entropy  ds  *  o  and  CpdT  s  T  lUl  dp 

aT  p 
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i  divide  second  into  first 


CvdT  =  -T  dv 

3T  V 


(^)  (_^)  (3£) 

3T  p  3p  V  3v  S 


but  from  the  mathematical 


theorems  we  may  rearrange  our  second  term  of 


(3p)  (3v)  .  ( 3p) 

TvT  3Tp  3Tv~°  equality  for 


(3v)  (3T) 

3T  p  3p  V 


(3P>  (3v)  .  Op) 

3v  T  sT  p  7t  V  * 


(^)  (3£)  (3V)  _  _  ^  (3V)  (3T) 

ap  V  av  T  3T  p  “  ”  *  3T  p  3p  V 


( 3v) 

3P  T  * 


Substituting  in  the  specific  heat 


ratio 


C£ 

Cv 


equation  we  get: 


.  ,  (3T)  (3£) 

Cv  aT  p  aP  v  av  s 


(^)  (3^) 

ap  T  av  s 


Ov)  ( 3p) 
ap  T  av  s 


or 


Cp  _  ( 9v) 
Cv  ■  ap  T 


(ap) 
av  s 


There  is  possibly  another  relationship  between 
P,  V,  and  T  for  entropy  so  that  the  fundamental  Tds  equation 

might  be  used  for  ,  where  an  identity  might  be  obtained 
t  X  Cv 

for  ,  This  may  be  accomplished  by  using  a  fundamental 

av  s 

mathematical  theorem  which  will  relate  entropy  (s)  to 
pressure  (P),  volume  (V),  and  temperature  (T).  It  is  known 
that  entropy  is  a  function  of  p,  v,  and  T  and  that  a 
relationship  exists  among  p,  v,  and  T  so  that  entropy  may 
be  regarded  as  a  function  of  any  two  of  p,  v,  and  T, 
Similarly,  any  one  of  p,  v,  and  T  may  be  considered  a 
function  of  entropy  and  one  other  of  p,  v,  and  T.  Thus, 

regarding  p  to  be  a  function  of  s  and  v, 

dp  s  (3,^^  ds  ♦  ( ^P^  dv  ’,  and  considering  v  to  be  a 
3s  V  av  3 
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function  of  s  and  T,  dv  =  ds  +  dT  ;  substituting 

3s  T  3T  s 


this  expression  for  dv  into  the  former  equation: 
dp  =  ds  ♦  <|£)^  ds  *  <|i)^  dl] 


ds  +  dT 


3T  s 


dp  s  ds  +  ds  +  dT 

3s  V  3v  s  ds  T  3v  s  3T  s 


dp  =  (il)  ]  ds  +  (^)  (3v)  dT  . 

^  3s  V  3V  s  3S  T  3V  s  3T  s 


Equating  the  dT  terms  of  the  equations  for  dp  we  get: 


“P  -  T  *  It 
dp  =  r<i£>  ♦  (SE> 


r<i£>  ♦ 

(  3P) 

t1  ds 

.  r(i£) 

(3v)  ] 

1  JS  V 

3V  s 

3S  Tj 

1  3v 

s  TT  sj 

.  (|E) 

3v  S 

(  3v) 
3T  s 

or 

1  .  < 

3T)  (3£) 

3p  S  3V 

(i£) 

(^T) 

so 

St. . 

(^>  (^) 

*  3T  s 

av  s 

Cv 

3p  T  ^ 

s  3V  s 


From  the  equation  of  the  difference  in  specific  heats, 
we  will  obtain  the  relationships  of  p,  v,  T  and  other  con¬ 
stants  for  three  types  of  gases.  The  first  will  be  the 
ideal  gas  equation,  pv  :  RT 

(RT)  ^ 

Cp-Cv  ^p-  3  T  p  *  p2 


‘I^’t  •  -  ”2  *  f  1!  -It  ■■ 


.  R iTi 


e2  'p2 

=  1 


P'  P*  V' 


Cp-Cv  =  R  for  an  ideal  gas. 
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Next  we  will  get  the  difference  in  specific  heats  for 
a  gas  that  follows  the  Van  der  Waals  equation  of  state. 

(P  +  a)  (v-b)  =  RT  ;  T  =  i  (p  +  a)  (v-b) 


v2 


V' 


fv’p  =  ?  [*'’  *vT  *  ‘^’1  !  squaring  this  term 

i-  f<P  >  '  (''->>>  f2*>l  ^  *>,• 

e  get:  ^2  [  ^2  “3J  J  inverting  this 


we  obtain; 

aT  p 


C3£)  _  2a _ RT 

av  T  " 

Cp-Cv 
Cp-Cv 


R2 


[« ■  ^1 

2a  (v-b) 2  -  rtv3 


-  2a(v-b)^ 


<v-b)2 

(v-b)2 

L  (v-b)2 

T  <3v)2 
^  IT  P 

ISv  T 

TR2 

1 

RTv*  -  2a  (v-b) 2 

fP  ♦  «  - 

1  y2 

2a(v-b)l2 

'"  vr- J 

v^  {v-bV=^ 

3 


«  R*T 


1 

rRT  v^  -  2a(v-b)* 

fP  -  2a  (v-b)|* 

1  v*(v-b)* 

1  V2  yi  J 

1 

Van  der  Waals'  equation* 


The  difference  in  specific  heat  for  a  gas  that  obeys 
the  Callendar  equation  of  state: 


RT  s  p(v-b  ♦  a  )  where  n  «  in 
jn  » 
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■  (RT  . 

(^)  _  a _ £_  -  +  b) 

ST  p  '  3  T 


"  P  +  1 


(  3v)  2 
aT  p 


R  ^  na _  (^) 

p  +  1  *  av  T 


(v-b  +  Tr>)  0 
(v~b  ■♦•  a)^ 
Tn 


cp-cv  =  -T  =  T  f*  *  2a 

^  3T  p  av  T  P  Tn 


P  Tn  +  IJ  I  (v-b  + 

*■  -j-n 


=  RT2 


RTn+1 


+  pna 


Tn+1  <v-b  ♦  a  ) 
**  <j>n 


^  V  10  .  r,  4.  1  -  ii. 

;  where  n  = 

3  ^ 


13  10 

Cp-Cv  S  RT^  £2X— £IX-S_ 

pT  (v-b  ♦a  ) 


Callendar's  equation. 


Berthelot's  equation  of  state: 


p  =  £I  -2-, 
v-b  Tv‘ 


RT _ a 

v-b  Tv 


-:Zz  -  P  -  ^ 


Cp-Cv  s 

3T  p  3v  T 


av/T 


[  (RT  -a  ) 
.  1 3  v-b  Tv^ 


-RT _  ^ 

(v-b)2  Tv^ 


(RT  -a  -  p) 

f.(RT) 

av^  Tv^  - 

*  ^ 

MTv^/ 

-  r3£] 

aT  p 

* 

aT  p 

s 

aT  p 

laTj 
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0  =  (v-b)R  -  RT  (dv) 

_ ar  p 

( v-b)  2 


r  (dv)  1 

(Tv2)  0  -  a  It^v  ar  p  +  v^J  -  0  =  0 


(Tv2)2 


0  =  R  (v-b)  -  RT(^) 
_ dT  p 

(v-b) 2 


2aT  v(dv)  +  av2 

- . £■..,_ . 

(Tv2)2 


0  =  RT2v'*(v-b)  - 


RT3v'*(^) 
dT  p 


2aTv(v-b)2(dv)  +  av2(v-b)^ 
dT  p 


TV*(v-bK 


RT2v'*(v-b)-RT3v‘*(^)  +  2aTv(v-b)2  (^)  +  av^(v-b)^  =  0 

dT  p  dT  p 


RT^v‘*(v-b)  +  av^  (v-b)^  =  (^)  RT^v**  -  2aTv(v-b)2 

dT  p  I 

(^)  .  RT^v**(v-b)  •»•  av^ 

3T  p  r'tM  2aTv(v-b)^ 


RT^v^(v-b)  av  (v-b)2 
RT^v^  2aT(v-b)^ 


(  3v)^ 
3T  p 


RT2v‘*(v-b)  ♦  av2(v-b)2 

RT^v^  -  2aTv(v-b)^ 

»  , 

Cp-Cv 


Cp-Cv 


^  (dv)2  (^) 

3T  p  T 


RT2v3(v_b)  ♦  av  (v-b)2 

2 

“£i  ♦  is. 

rt3v3  -  2aT(v-b)2 

k 

(v-b)2  Tv3 

11 


12 


a 

RTF  +  1 

v-b  [rW 

+  l] 

T(v-b) 

(v-b)^ 

[  (v  -b) 

T 

or 

a 

RTv  +  1 

(^)2 

a 

1  +  RTv 

- 

r  a  '  _ 

3T  p  ‘ 

r  a 

JLl 

T  [rTv2  ” 

v-b  J 

T  LRTv2 

v-bj 

Cp-Cv  = 


[  — 1 

-  .  1  •»•  RTvJ 


T  [rWT  “  v-b] 


-a 

RTv 


RTa 


[l  ♦  RTv  j 

2 

1  a 

r  -  a  1  • 

LrTv2  *  v-b. 

2 

<v-b)2  RTv2  (v-b) 

RTv 2 (v-b) 


-a 

RTv 


Re 


_ 1_ 

(v-b) 2 


Using  properties  of  superheated  steam  from  a  point 
near  the  saturated  line  to  a  point  far  away  from  the 
saturated  line,  we  may  compare  the  values  of  these  three 
expressions  for  the  difference  in  specific  heats. 

The  first  to  be  calculated  will  be  the  id^al  gas  case. 
This  is  a  constant  and  is  equal  to  the  universal  gas  con¬ 
stant.  The  universal  gas  constant  is  approximately  1545 
ft  Ib/lb  mol«R.  Cp-Cv  *  1545  ft  Ib/lb  mol®R  for  the  ideal 
gas. 


To  solve  for  the  two  constants  we  will  employ  the 
properties  of  the  critical  isotherm  at  the  intersection  of 
the  saturation  line  for  the  Van  der  Waals,  Berthelot,  and 
Dieterici  equations.  For  the  callendar  equation,  we  will 
choose  two  points  in  the  superheated  steam  range  of  values 
to  obtain  these  constants.  This  is  done  because  the 
Callendar  equation  does  not  apply  at  the  critical  point. 

Berthelot*8  equation; 

Since  the  critical  point  is  the  limiting  position  as 
t%io  points  on  each  of  the  saturated  liquid  and  the  saturated 
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vapor  at  the  same  pressure  and  temperature  approach  each 
other,  it  follows  that  at  a  maximum  point  on  the  curve,  the 

slope  is  zero.  This  means  that  =  Oj  also,  since  the 

critical  point  is  a  point  of  inflection  on  the  critical 

isotherm,  ^y-^2T  ' 

From  these  two  facts,  plus  the  equations  themselves,  we  may 
determine  the  two  constants. 


'-2a  \  -  2a  - 

(v-b)2 

av^T  ~  Tv**  ^  (v-b)  ^  Tv^  (v-b)  2 

6a  2RT  2a  _  RT  ^  .  RT 

TviT  =  (v-b)3  »  Tv3  ■  (7Ib)2  Tv**  ■  (v-b)3 


Berthelot’s  equation;  p  =  ^2 


(^)  /av^  \  -  (aT^A  ,  -  RT 

av  T  ■  \  a  V  T/  \  av  /T  (v- 


O^p)  -  6a  ^  2RT 


”2 


2a  RT 


Dividing  one  by  the  other,  we  eliminate  a  and  obtain  b, 
in  terms  of  the  critical  volume. 


RT(v-b)^  . 
Tv 3  3a  (v-b)2  RT 


2 

7  V 


2 

v-7v  =  b;  or  v^,  =  b 


3 

Putting  this  into  the  first  equation  we  get: 

2a _  _  RT _ 

T(3b)3  ‘  (3b-b)2  * 
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RT 

(v-b)2 


»a 


-a 


e  RTv 
(v-b) 


0 


15 


16 


0 


1 _ L-  .  f-i-  -  1  ]\ 

tv-b) 2  (v-b)'^  '  iv-b  v-bj 

1-2  .  (v-b)''  .  V 

Tvrn’2  "  *  “TTTF  =  2  = 


2v-2b  =  V  }  2v-v  s  2b  ;  v  s  2b  . 


Substituting  this  back  into  the  first  equation  for  derivative 

we  get:  -  -j;.  ..  =  0 

RT(2b)2  2b-b 


a  1  ab  a 

R^'Ubz  =  F  Ri+bz  "  Tc  = 


Substituting  these  answers  back  into  the  equation  of  state 
we  find: 


Pc 


Ra  e 


a4Rb 

l?a5F 


4Rb(2b-b) 


Van  der  Waals*  equation  of  state: 


P  s 

RT _ a 

V-b  ”  ^ 

/3P\  .  -RT  2a  .  - 

\7v/  T  ”  (v-b)^ 

2a 

RT 

v^ 

*  (v-b)2 

-  2RT  .  6a  .  g 

3v^  T  ”  (v-b)^  ^ 

3a 

vJT 

RT 

“  (v-b)  3 
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2av‘*  _  RT  (v-b)3  _  2 
v^Sa  '  ’(v-b)^  RT  "  3 


2  1  2a  RT  2a  RT 

V  -  -  V  =  b  =  j  VC  }  -^^3  =  (3b.b)2  =  27ET  *  irb2 

2a  RT  _  8^ 

27b  ■  4  i  =  27Rb 

Finally,  substituting  into  the  equation  of  state  we  get: 

.  8a  .  _ a  _  4a-3a  _  a  _ 

27b(2b)  9b^  27b^  Tt^  "  ° 


Since  the  Callendar  equation  of  state  does  not  apply 
at  the  critical  point,  we  may  obtain  the  approximate  values 
of  its  constants  by  substituting  values  of  P,  T,  and  v  in 
the  superheated  steam  range  where  this  equation  is  applica¬ 
ble. 


Point  one  will  be  a  pressure  of  40  psia,  a  temperature 
of  300®F,  and  v  =  11.04  ft^/lb  specific  volume.  The  first 
point  will  be  p  s  80  psia;  T  s  500®F;  and  v  -  7.02  ft^/lb. 


as  (-V  ♦  b  ♦  RT)  T"  s(-l26.2  ♦  b  ♦  1545x960) (. 96 )  ^  x  10^° 
p  11,520 

a  =  873x10^  (b  ♦  128.5  -  126.2)  =  873x10^  (b  +  2.3) 

b  s  873x10^  (b  ♦  2.3)  -  4.8  s  8J73  b  ♦  873(2.3)  -  4.8 
401x10^  401  ‘♦01 

2.175b-b  s  4,8  -  5  s  -  0.2  =  1.175b 

b  s  -.17ft3/lbmol 
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Using  another  point  of  p,  v,  and  T  we  find: 


a,b,  p  =  160  psia,  T  =  900®F,  v  =  5,015  ft^/lb 

a  =  (-V  +  b  +  =  (-90  ♦  b  +  1545x136, 0)  (1, 36 )"^xl0^  ° 

P  23,050 

a  =  (b  ♦  91,1  -  90)  2,781x10^°  =  (b  +  1,1)  2,781x10^° 

Equating  these  two  we  get:  (b  1,1)  2,781x10^°  = 

(b  +  2,3),873xl0^°;(b  +  1.1)  2,781  =  b  +  2,3  ;  2,781 ^  ^ 

,873  ,873 

1,1(2,791)  =  b  +  2.3 

0,873 

3.185b  -  b  =  2.3  -  3.5 

2.185b  =  -1,2  b  =  -0,55  Callendar’s  equation. 


a  =  (b  +  1.1)  2.781x1010  =  (la  -  0.55)  2. 781x101°  =  ia3xlolO 


For  the  Berthelot  equation  of  state  for  steam: 

Vq  =  0,906  ftVlbmol  Pc  *  459,000  Ib/ft^ 

Tc  =  705OF  =  1165®R 

b  =  Vc  =  0.906  =  0.302  ftVlbmol 
T-  3 


(12b  pc)2  s  2^  ; 

3b 

a  s  (.906)  278x10*” 
3.09  X  10 5 


a  s  3b(12bp  )^ 
2R  ® 


(.906)(4(. 906)4. 59)^x101° 
2  x  1545 


s  81.4  X  10^  =  8a4  X  10®  lb _ 

ft2(ft'’)2  OR 

Ibmol 
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a  = 


8.14  X  10® 

Ibmol 


OR  lb 
ft  2 


Calculation  of  Cp-Cv  for  Berthelot  equation  at  point  one, 


p  =  1440  psfa,  T  5  660°R 

V  s  700  ft®/lbmol 

a  =  8.14  X  10®  b 

(Ibmol) 2  ft 2 

Berthelot *8  equation; 


0.302 

Ibmol 


R 

v-b 

f2^\ 

Cp-Cv  =  T  - 

J  - 

RT 

(v-b)* 

- 

Tv® 

6.6  X  10* 

1545 

,  8.14  X  10*  1 

.700-. 

302  (6.6)2x492x10“] 

1545  X  660 

.  2(8 

.14)(10«) 

(700-. 302)2 

(660)343x10® 

Cp-Cv  •  6$0  fa.»0n»3  ♦  0.00381p 

2.081  -  0.00719 


Cp-Cv  =  1556 

Dieterici*8  equation  of  state; 
p  =  1440  psfa  \  T  s  660®R,  v 

a  s  3.2x10®  ;  b  s  0.453 


700  ft® 
Ibmol 


RTvf  -  * 

Cp-Cv  s  Ra  [1  ♦  RTv 

—  1* 
RTv*  v-b  J 


_ 1_^  .  a 

(v-b)*  RTv* (v-b) 
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-3.2  X  10^ 


Cp-Cv  =  1.545  X  10^  e  1.545  x  10^  x  .66  x  10^  x  7,  x  10^ 


3.2  X  10® 


1.545  X  10^  X  .66  X  10^  x  49  x  10** 
3.2  X  10® 


700 


■  3.2 

1  +  1.545  X 


■■  2 

X  7  X  102 


[49  X 


10' 


103  X  .66  X  103 

3.2  X  10®  ^ 

1.545  X  103  X  .66  X  10^  x  49  x  lO**  x  7  x  lO^J 


Cp-Cv  = 


1.545  X  103  [1  +  ,004] 


.0045  r _ 

e  L15.6 


X  10“ 


- 

7  X  102j 


r  1  _ _ 

[49x10“  109 


i - 1 

9.xlO®J 


Cp-Cv  s  1.545  X  10^  Cl. on  [.0204  x  10 


.  0092  X  10~^  ] 


[1.01-]  [.641  X  10 


—  4 


,143  X  10-2]2 


1.545  X  10^  2.03  X  10-* 

(-0.1374  X  10-2)2 


1545  [2,03]  X  10-® 

(1,37)2  y  10“® 


1545(2,03) 

1,88 


1668  ftlb 


lbmol®R 


a  =  1400  atm  (ft^  )2  ;  b  =  0.488  ftVlbmol 

Van  der  V/aals  equation. 


1400  X  2116.2  ^  (ft^  )2 

ft 2  Ibmol 


2,96  X  10®  (ft^  )2 

ft  2  Ibmol 


T  s  200®F  =  660®R;  v  *  38.85  =  700  ft^ 

Ibmol 
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R  =  1545  ftlb 


p  =  10  psia  =  1440  abs, 

lbmol°R 

Cp-Cv  for  the  Van  der  VJaals  at  this  point  will  be  calculated 
first. 


Cp-Cv 


1 — 

CM 

2 

‘rTv3  -  2a(v-b)2' 

P  +  -  2a(v-b) 

.  V  2  V*  -J 

(v-b)2 

s  (1.545)2x10^ 


660 


1440  +  2.96  X  10^  -  2(2.96  x  10^(7.  x  102) 


L 

(7. 

)2  X  lO** 

73 

X  10®  J 

X  [1.545  X  10 

^  X  6.6 

X  10^  X 

7^  X  10^ 

-  2(2.96  X  10^)72  X  lO*"! 

L  73 

<  10®  X 

72  X  10 

4 

j 

=  2.39  X  10® 

660 

1/2 

- 

2 

1440  + 

296  -  6 

X  10® 

m 

49 

49  X  lO**. 

350.0  X  10>2 

-  300 

X  10‘0 

1  n  6 

660  V2  2 

17,500  X 

10‘0 

~  j  3  X 

xu 

1440  +  6.05  -  12. 2_ 

[3471  [660  1  [3471 

LlnJ'  [rTTj 


347lx  10®  =  2.18  X  10®  ftlb 

1.436  X  103  lbnol°R 


-  ftlb 

lbnol°R 

Cp-Cv  for  the  Callendar  equation  for  superheated  steam  at 
these  same  conditions  is: 


22 


2 


Cp-Cv 


RT2 


1  3 


pT 


3  pa 

(v-b  +  a  ) 


Cp-Cv  =  1545(660)2 

1  3 

1545(660)  ^ 

10 

+  “T  (1440)(2.96x106) 

1  0 

1440(660)  ^ 

(700-0.488  +  2.96x10® 

1  0 

(660) 

Cp-Cv 


1.545x(6. 6)2x10"^ 


I  3  1  3 

1.54SxlO^(.66)  ^  x(103)  ^  ■>•1.44x10^° 

II  11 

1.440xl0^(  .  66)  ^xdO^)  ^  (7x102»2. 96x10^ 

1  0  1  0 

(.66)  ^x(103)  ^ 


Cp-Cv 


1.545x(6.6)2  X  10^ 


0.2558  X  10^^ 
0.2385  X  10^® 


12 


1.44  X  10^° 

(7  X  10^  +  18  X  10 


,  4  ro.2558  X  10^^  Y 

Cp-Cv  =  1545  X  (6.6)^  X  10  Li.6506  x  10'“  J 

Cp-Cv  s  1545  X  43.6  x  10*’  [0.1545]^  x  10“** 
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Cp-Cv  =  1545  X  43.6  x  0.0239  =  1545  x  1.04 


=  1610  ftlb 

lbmol°R 


Cp-Cv  =  R  =  1545  ftlb  =  constant  for  ideal  ^as. 

lbmol°R 

To  see  how  closely  these  gas  equations  represent 
superheated  steam,  we  will  take  two  points  for  each  gas  to 
determine  the  change  in  entropy  for  a  reversible,  isothermal 
expansion  of  steam. 

Fiom  the  equation  for  the  change  in  entropy  as  a 
function  of  two  thermodynamic  variables,  we  get: 


ds  = 

(3s) 

dT  + 

(as) 

dv  and 

II 

3T  V 

3v  T 

3V  T  8T  V 

ds  = 

(as) 

dT 

( ap) 

dv.  For 

a  reversible  isothermal: 

3T  V 

3T  V 

dT  = 

0  so 

ds  = 

(3^,) 

dv. 

3T  V 


(^) 

3T  V 


for  a  Van  der  Haals  equation: 


P 


(RT 
a  v-b 


3T 


R _ 

v-b 


Callendar  equation;  p  = 


RT 


v-b  ♦  a 


24 


(RT _ ) 

v-b  +  a 
9  Tn 


3T 


Jv 


(v-b  +__a  )  R  -  (RT)(-na  ) 

_ 1 

(v-b  +  ^ 

ipn 


-a 


^ .  .  .  .  RT  RTv 

Dieterici  equation:  p  =  — i.  e 

v-b 


alRL-S. 


-a 

RTv 


v-b 


3  T 


-a 

RTv 


-a 

RTv 


-a 

RTv 


V 


[■  RTv  RTVI  RTv 

Te  (  a  )  e  .Re  r  a 
RTv^'*’  -I  v-b  [rvT 


perfect  f;as:  pv  =  RT 


p  .  £I 


3  (W 

V  ]  _  R 
3T  /v  V 


Bert helot  equation: 


_  ^ _ ^ 

v-b  Tv^ 


3(RT _ 

v-b  Rv^ 
3  T 


V 


R 

v-b 


♦ 


a 

T2v2 


ds 


dv  Van  der  \7aals  equation: 

v-b 


+  1 
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ds  = 


a  RTna\ 

(v-b  *  T’^)  R  \ dv 


(v-b  +  a  )  ^ 
fTT 


Dieterici 


7 


■  +  RT  na  dv 


(v-b  +  +  l(v-b  +  a  )2 

■pri  rpM 

Callendar’s  equation. 


-a 


ds  =  Re.Rtv  r  a  ^  dv  = 

v-b  [RvT  J 


D 

ds  =  —  dv  perfect  p,as 

V 


a  dv 


R  dv 


RTv 

e  (v-b)vT 


a 

RTv 

e  (v-b) 


ds  5  ^-£1  dv  +  -^ 
v-b 


T2v2 


dv  Derthelot*s  equation. 


To  get  the  change  in  entropy,  we  will  integrate  the 
equations  for  -  Sj, 


S2 

ds  5 


/•  2  P  1  2  r 

I  dv  =  R  In(v-b)  =  R  ln(v 

V  J  v-b  L  Jvj  L 

r  (V2-b)1 

•  ^r"(Vi-b)J 


,-b)-ln(  v, 


i-b)] 


Van  der  VJaals 


®2 

f  ds  =  R  f,— .121  ■-  -»•  J* 

Y  vy  v-b+a  T>^‘*’1  yf 


v-b+a 


dv 


(v-b^a  \ 

rnj 


2  - 


Rrin(v-bt|^)l'''  + 

*■  ■^''l 
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d  = 


2 

ds 


-a  -a 

RTv  RTv 

Re  a  dv  Re 

+  — —  dv 
(v-b)  RTv  v-b 


a 

‘  RTx 
a 

let  RTv  ■  ^ 


dx 


-a 

RTv  2 


dv 


dv  =  -RTv^  (jx 
a 


ae  RTv 2  ^  Re  C»RTv^  dx) 

Tav(v-b)  ^  *  (v-b)  a 


Re 


-X 

vdx 

(v-b) 


R^T  e~^  vdx 
a 


-X 

Re  dx 


R2Te~^a  dx 
a  RTx 


Re~^dx 

X 


let  v-b  ■  V 
where  v>>b 
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Perfect  gas* 

V2 

In  (v-b)  + 

[•  1 

L  Jv, 

> 

1 

s  R 


(V2"b) 

(vi-b) 


a  a 

.j,2y2  .p2^i  Berthelot’s  equation. 


The  two  points  used  for  the  test  for  four  of  these  gases 
will  show  how  closely  each  represents  the  change  in  entropy 
of  superheated  steam*  The  first  will  be: 

temperature^  =  ^00°F  =  1060°R;  pressure^  =  1200  psia  = 

173,0001b  .  =  0,4016  =  T.Sft^ 

ft 2  *  Ibmol 


The  second  point  will  be:  temp  - 
pressure  2  =  20  psia  =  28801b/ft2 

Van  der  V/aals'  equation:  s.  -  s 

■  '  I  ■"  a  I  ■  I  i-i  2 


600°F  =  1060°R; 

;  V2  =  31.47  =  626ft® /Ibmol 

(626  -  0*488) 

“  ^  Xil  p  A  I.  A  a 
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S-  -  s  =  1545  In  =  1545  x  In  (89.05) 

2  ^  7.012 

s  -  S  =  1545  X  4.485  =  6930  ftlb 
^  ‘  Ibmol 


Callendar's  equation; 


Sj  -  Sj  =  R 


.  T„. 


Vi-  b  +  ^  1 


a 


v,-b  +  —  v,-b  +  — 

1  <pn  ^  'p’™  j 


where  a  =  1.53x10^°  b  =  -0.55 


S2  -  5  R 


m 


(626  ♦  .55  +  1.53x10^°  \ 

_ (1.06x103)H  )+  1060(10)  1.53x10^° 

7.5  ♦  .55  ♦  1,53x10^°.  .  / 

.  i  g 


(1.06xl03)-y 


3  (1. 06x10  3  )-*•! 


7.5  ♦  0.55  +  1.53x10^°  626  +  .55  +  1.53x10^° 

(1.06)-410^0 


(1.06-^)10loJ 


S2  -  Si  =  R 


,  (626.55 

♦  i.26)  .  . .  r 

1 

1  ll 

8.05  ♦ 

1.26  *  "•""L 

9.31 

627.8lj| 

Sj  -  Sj  a  R  |4.21  ♦  0.476}  a  1545  }4.686| 


S,  -  Si  a  7240  ftlb 
^  TOnol 
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a  3.2x106 

RTV,  °  1.545x103x1.06x103(7.5)  =  Dleterlci's  ..juatlon. 


£ _  ^  0,26x7.5 


RTv^  V2 


626 


0.0031 


s,  -  s,  =  1545  .0031  - 


1 


[• 


.  00  3  1 


-  In  (.0031)  -  0.26 


*  ^^6  +  ln(.26)  +  .015] 


[• 


Sj  -  Sj  =  1545  |.0031  -  1  +  5.76  -  .26  +  1.3  -  1.345  +  .OI5] 


S2  -  Sj  =  1545  X  4.47  =  6900  ftlb 

Ibmol 


Perfect  gas; 


Sj  -  Sj  s  R  In 

Vl 


1545  X  In  s  1545  x  4.42 
7.5 


®2  "  ®l  ■  6830  ftlb 
Ibmol 


Berthelot's  equation; 


b  s  0.302;  a  =  8.14x10®  lb  (ft®  )  oj^ 

Ibmol 


s,  -  s, 


R  in  ,  a 

v,-b  5*77 


S2  -  Si  s  1545  In 


(626  -  .302)  _  8.14x10® 


7.5  -  .302  (1.06)2x106x6.26x102 
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^  8.14x109 _ 

(1.06)2x109x7.5 

Sj  -  Sj  =  6900  -  1.158  +  0.967  x  102 

=  6900  +  95.5  =  6995.5  ftlb 

Ibmol 


From  the  steam  tables  by  Keenan  and  Keys,  "Thermodynamic 
properties  of  Steam,"  we  obtain  the  actual  change  in  entropy 
of  the  two  states.  Entropy  at  20  psia  and  600®F  =  1.940, 
and  entropy  at  1200  psia  and  600°F  =  1.410  Btu/lb. 

This  is  equal  to  1.940  -  1.410  =  0,53Btu/lb. 


In  ftlb  ,  this  is  0.53  x  18  x  778.3  =  7410  ftlb  . 
Ibmol  TEinol 


The  variation  for  each  gas  equation  is  shown  in  Figures  1 
and  2 . 

The  percentage  deviation  for  each  is  as  follows; actual-calc. 

actual 


Van  der  Waals  -  7410-6930  _  480  x  100%  =  6.48% 

7410  7410 

Callendar  -  7410-7240  s  2.29% 

7410 

Dieterici  -  7410-6900  s  6.87% 

7410 

Perfect  gas  -  7410-6830  =  7.83% 

7410 


Bert he lot  -  7410-6996  =  5.60% 
7410 
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It  is  obvious  that,  of  these  cases  considered,  the 
Callendar  equation  describes  superheated  steam  best.  The 
same  could  be  done  for  other  gasses  to  determine  which 
equation  would  fit  them  best. 

Ill,  CONCLUSIONS 

It  can  be  seen  that  the  differences  and  ratios  of 
specific  heats  and  change  in  entropy  may  be  determined  for 
the  considered  five  equations  of  state. 

Comparisons  of  the  difference  in  specific  heats  and 
change  in  entropy  for  superheated  steam  at  a  few  selected 
points  show  that  these  equations  fit  the  actual  case  fairly 
well,  and  most  of  the  equations  are  useable  for  engineering 
calculations. 

These  reduced  equations  could  be  used  to  establish  use¬ 
ful  engineering  design  tables  for  several  frequently  consid¬ 
ered  gases. 
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Figure  1.  CHANGE  IN  ENTROPY  OF  SUPERHEATED  STEAM  CALCULATED  FOR  FIVE  GAS  EQUATIONS. 

T=G00OF  CONSTANT 


aSuvqo  Xdoa^us  03  pasn  uo'(:)*nb3 
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Figure  2.  PERCENT  DEVIATION  FROM  ACTUAL  CHANGE  IN  ENTROPY  FOR  SUPERHEATED  STEAM 
FOR  FIVE  GAS  EQUATIONS.  T=600OF. 
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